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Charmonium in lattice QCD and the non-relativistic quark-model 



OO 
O 

o 

(N 

Or 
<D 

in 



43 

Oh' 
i ' 

43 



> 

(N 
00 

in 

(N 

On 
O 
00 

O 

> 

X 



Jozef J. Dudek 1 ' 2 and Ermal Rrapaj 2 

(for the Hadron Spectrum Collaboration) 

1 Jefferson Laboratory MS 12H2, 12000 Jefferson Avenue, Newport News, VA 23606, USjQ 
2 Department of Physics, Old Dominion University, Norfolk, VA 23529, USA 

We compare the results of a numerical lattice QCD calculation of the charmonium spectrum 
with the structure of a general non-relativistic potential model. To achieve this we form the non- 
relativistic reduction of derivative-based fermion bilinear interpolating fields used in lattice QCD 
calculations and compute their overlap with cc meson states at rest constructed in the non-relativistic 
quark model, providing a bound-state model interpretation for the lattice data. Essential gluonic 
components in the bound-states, usually called hybrids, are identified by considering interpolating 
fields that involve the gluonic field-strength tensor and which have zero overlap onto simple cc model 
states. 
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I. INTRODUCTION 

New results from CLEO and the B-factories have mo- 
tivated much new interest in the spectrum of charmo- 
nium. Several new resonances above open-charm thresh- 
old do not appear to fit in the patterns established by 
the clear sub-threshold states. These lighter states ap- 
pear to be well described as cc states bound in a po- 
tential which confines as large distances, but within this 
framework it has proven difficult to explain many of the 
new resonances leading to suggestions that we should in- 
clude physics beyond that allowed in the simplest po- 
tential models, specifically that we require a gluonic 
field carrying non- vacuum quantum numbers ( "hybrids" ) 
or higher Fock states such as cuuc ( "multiquarks" or 
"meson-meson molecules" ). The clearest signal of physics 
beyond the quark-model, known as an "exotic" , is when 
a state has external quantum numbers not accessible to 
a cc pair, be that in the J PC , such as 1 ^ or in the 
flavour, such as an isospin 1 state with strong decays to 
charmed mesons 0, Q ■ 

Much of the phenomenology of the old and new char- 
monium is performed with quark potential models and 
their extensions, such as the Coulomb gauge model [||, 
or the flux-tube model These models usually work 
within a Fock state basis, fixing the number of con- 
stituents which have a prescribed set of interactions, usu- 
ally motivated by, but not formally derived from QCD. 
The models often predict a mass spectrum and possibly 
radiative and decay properties which can be compared 
to experimental data. Assignment of experimental state 
structure follows if the theory matches the data. 

An example of this kind of phenomenology would be 
the vector (1 ) channel of charmonium in which es- 
tablished resonances at 3097, 3686, 3770 MeV and en- 
hancements at 4040,4160,4420 MeV are identified as 



the model cc states with internal quantum numbers 
n 2S+1 Lj = 1 3 S 1 ,2 3 S 1 ,1 3 D 1 ,3 3 S 1 ,2 3 D 1 ,4 3 S 1 0. The 
subsequent suggestion of another vector state at a mass 
4260 MeV[l, 0,H, Q is thus supernumerary within this 
scheme and proposals have been made that one should ac- 
count for it by considering degrees-of-freedom not present 
in the quark potential model, such as explicit gluonic 
quantum numbers, making this a hybrid state, or pos- 
sibly the quantum admixture of a conventional and a 
hybrid state Since the state has external quantum 
numbers accessible to a cc we choose not to label it ex- 
otic, with "crypto-exotic" implying a hidden exoticness. 

Distinguishing between "conventional" and "crypto- 
exotic" states is inherently difficult within lattice QCD 
calculations; unlike in bound-state models where one 
must specify upfront the structure of the state, in lat- 
tice QCD one simply provides interpolating fields, con- 
structed from quark and gluon fields, of the appropri- 
ate external quantum numbers and extracts mass infor- 
mation from QCD correlators. Structure information is 
embedded in the relative magnitudes of overlap of var- 
ious interpolating fields on to each state, but it is not 
always clear how one should interpret this. In this pa- 
per we propose a model-dependent scheme based upon 
a non-relativistic bound-state problem motivated by the 
physics of quark potential models commonly used to de- 
scribe heavy quarkonium. 

In [TlT ] , a set of interpolating fields 1 built from fermion 
bilinears with up to two discretised covariant deriva- 
tives was used to extract a meson spectrum with charm- 
mass quarks. The spectral quantities extracted where 
the masses of the states (ttiat) and the overlap of 
each state onto the particular interpolating fields used 
(Z[ N) = (0|e>i(0,0)|iV)) as expressed in the spectral de- 
composition of a Euclidean two-point correlator at zero- 
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we will also refer to these as "operators" 
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where the sum over states labelled by N extends over all 
meson states with the quantum numbers of the interpo- 
lating fields Oij. Details of the method used to extract 
771, Z can be found in (Tlj . 

In this paper we compute these overlap factors us- 
ing a non-relativistic reduction of the continuum limit 
of the interpolating fields and cc states constructed as in 
a quark potential model. Comparison with the numeri- 
cal values extracted from the lattice calculation tests the 
potential model formalism and gives a model-dependent 
interpretation of the lattice spectrum that agrees reason- 



ably well with the conventional quark-model picture, but 
goes beyond it by suggesting the influence of gluonic de- 
grees of freedom. 



II. NON-RELATIVISTIC POTENTIAL MODEL 
STATES 

In the non-relativistic potential model we construct 
states of definite quark and anti-quark number and as- 
sume that the only effect of the gluonic field is in provid- 
ing a static potential which binds the quarks into mesons. 
A generic unflavourcd qq meson state can be constructed 
which is an eigenstate of orbital angular momentum L 
and total quark spin S, 



\n 2S+1 Lj,mj;d=y/2E^ £ (Lm L ; Sm s \ Jmj) ]T(ir; ±s\Sm s ) fA^^^Y^ (q) a\ q)b\ (±j?- $)|0), 

rriL,ms r,s * ' 

(i) 

where a\(p) / b\(p) is the creation operator for a quark/ antiquark of momentum p and z-component of spin r. The 
momentum-space wavefunction ^ni(|^1)j which carries the orbital angular momentum quantum number L and a 
principal quantum number n, is normally determined by solving a Schrodingcr equation with a potential of phe- 
nomenological origin . Within such models, states of different L can be mixed by non-central interactions such as a 
tensor force which are usually considered to be relativistic corrections to the dominant central potential, suppressed 
by powers of \q\/m q . 

Although this state is constructed to have non-zero total momentum, it does not transform in a Lorentz-covariant 
manner under boosts, the model having only Galilean invariance - we will see later that this reduces the usefulness 
of the model away from states at rest. For the bulk of this paper we will consider only meson states at rest and will 
leave the momentum-space wavefunctions unspecified. 



III. NON-RELATIVISTIC REDUCTION OF 
INTERPOLATING FIELDS 



A set of operators was presented and used in , based 
upon an extension of the set proposed in [I2| |. These 
operators used a simple discretisation, 

^jf( x ) = &(Uj(x)f(x + fa)-U](x-ja)f{x-ja)) 
-> D j f(x) + 0(a 2 ) 

of the covariant derivative and were constructed to trans- 
form irreducibly under the group of rotations allowed on 
a cubic lattice. At zero momentum there are five ir- 
reducible representations, A±, Ti, T2, E, Ai in which the 
various continuum spins are distributed (l3j . The opera- 
tors used were constructed such that although they trans- 
form irreducibly under lattice rotations, they also have a 
continuum limit in which they overlap with only a single 



state 2 of definite J PC . These operators then are expected 
to have an "unsuppressed" overlap with one particular 
J PC along with overlaps with "lattice artifact" states, 
suppressed by powers of the lattice spacing a. For exam- 
ple, the operator ijj\7 itp has a continuum limit tp D itp 
which overlaps with only f at zero momentum. At 
finite lattice spacing this operator transforms as T\ and 
hence can have overlaps with 3 ,4 ... suppressed by 
at least one power of a. 

We intend to compute overlaps of the type 

Z = (0|^(0)r^^ . . . ^(0)|t7 2S+1 j L j , mj; 0), 

in the limit that the internal momentum of the quarks in 
the meson is much smaller than the quark mass. We will 
use the free-field expansion of the quark field operators, 



2 in a few cases there are two continuum overlaps, see the appendix 
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and, since the quark-model state contains no gluonic field 
operators, we will neglect the gluonic field part of the 
covariant derivative 3 , Di = di + A, — > <9;. In many-body 
approaches to QCD, such as the Coulomb-gauge model 
[3], neglecting the transverse part of the gluonic field as 
a first approximation in this way is natural as the gluons 
appear with an effective mass such that their involvement 
in the spectrum is suppressed by a mass gap. 

In general the overlaps are found to be proportional to 
an object 

(2tfO (2i£)- v s (-q)T Ur (q). 

The spinor contraction is evaluated for the 16 possible 
gamma matrices using an explicit representation and ex- 
panded in powers of . The integrals over angle can be 
performed by expressing the components of q in terms 
fY™(q)'s, and the resulting expression, when combined 
with the quark spin and spin-orbit Clebsch-Gordans of 
Eqn[T]gives rise to overlap with only certain quark model 
states. The overlap is expressed via a single remain- 
ing integral over q = \q\ featuring the unspecified radial 
momentum-space wavefunction, <p n i,(q). 



A. Quark Smearing 

In [rH |. smeared quark fields were sometimes used - 
these have the continuum limit exp \_\<? 2 ^ D i D i \ip(x), 
where a is called the "smearing radius" . Considering 
as a first approximation the covariant derivatives to be 
conventional derivatives using these fields gives overlaps 
proportional to 

d*q 
(27T) 

x(2iq l ) (2iqj) 



v s (~cf)Tur{q). 



The purpose of smearing the quark fields is to enhance 
in a correlator the contribution of the ground state as 
compared to the contributions of excited states. Within 
the non-relativistic model we can easily demonstrate how 
this comes about - consider as an example the pseu- 
doscalar states which in this model have a n 1 ^ com- 
position. The interpolating field ip^f 5 ip has pscudoscalar 
quantum numbers and performing the non-relativistic 



reduction we find for the overlap with an n 1 ^ state, us- 
ing unsmeared quark fields 

Z(«> = ^2^/0 W (,)(l + ^ f ), 
and using smeared fields 



s(q) 1 + 



The unsmeared expression, if Fourier transformed into 
position space, is proportional to the value of the spatial 
wavefunction at the origin with a "relativistic" correction 
proportional to the second derivative of the wavefunction 
at the origin divided by the quark mass squared. To get 
a picture of how smearing might work to enhance the 
ground state over radially excited states, consider a sim- 
ple model in which we neglect the relativistic correction 
and use the bound state wavefunctions of the harmonic 
oscillator potential 5 . In such a model we can perform the 
above integrals exactly to find 
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so that if we smear with a = (3~ x we have overlap with 
only the ground state and with none of the excited spec- 
trum. The dependence upon a (3 for the lowest four states 
is shown in the inset of Figure [1] For charmonium a 
more realistic basis of wavefunctions are the solutions 
to the Schrodinger equation with the Cornell potential, 
— | ^ + fir - numerically solving this Schrodinger equation 
with the parameters given in [lij we find the overlaps 
shown by the curves in figure [1] It is clear that there 
is still a region of a in which the excited state contribu- 
tions are much suppressed relative to the ground state. 
Also shown are overlaps extracted from a lattice QCD 
calculation (with the same lattice parameters as those 
used in [111 ]) for various smearing radii 6 , where we see 
that there is generally reasonable agreement between the 
lattice smearing dependence and the model smearing de- 
pendence. 



IV. QUARK MODEL INTERPRETATION OF 
CHARMONIUM LATTICE DATA 

We have computed the overlap of the continuum, non- 
relativistic limit of the interpolating fields used in [ll[ 
with non-relativistic quark model states - the full results 
are tabulated in Tables HI HTllIIH These results should be 
supplemented with the fact that all of the B-typc opera- 
tors used in [ll| have zero overlap with any quark model 



3 we will discuss the consequences of this for operators involving 
the commutation of two covariant derivatives later. 

4 for the purpose of this argument we are assuming that the lattice 
is sufficiently fine that states do not resolve the cubic symmetry 
of the lattice and hence "lattice artifact" states can be neglected. 



5 We'll parameterise this such that the ground state ~ exp — A^-. 
for simplicity wc normalise such that the unsmeared ground state 
overlap is equal in the lattice and model calculations. 
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FIG. 1: Inset: Pseudoscalar overlaps with non-relativistic 
harmonic oscillator quark model states. Main Plot: Pseu- 
doscalar overlaps (including the "relativistic" correction) with 
numerical solutions to the Schrodinger equation with Cor- 
nell potential (solid curves) and lattice QCD simulation data 
(points). 



state, which follows from B ~ [D, D] — — > 0. A paren- 
thetical remark is that the operators with p and pi in 
their names can be added/subtracted to produce opera- 
tors with leading/subleading non-relativistic behaviour. 
This just corresponds to the usual notion of large up- 
per/small lower components of the non-relativistic limit 
of a Dirac spinor. 

Throughout this paper we use the following notation 
for the polarisation vectors of a spin-1 particle having 
momentum p and z-component of spin r: S^p, r). The 
polarisation tensors for spin-2 and spin-3 particles are 
similarly denoted. 

In what follows we will attempt a model-dependent 
description of parts of the lattice spectrum by compar- 
ing the relative size of overlaps extracted from a numer- 
ical calculation with the quark model expressions in Ta- 
bles U [TO [En] The lattice data used is an extension of 
that presented in pd| , using this time a somewhat larger 
basis of operators, in particular more operators of high 
mass-dimension. Doing this allows overlap with discreti- 
sation artifact states with a suppression that is smaller 
in powers of lattice spacing, a, so we would expect to 
see a somewhat more dense spectrum 7 . We remind the 
reader that we expect these lattice results to be domi- 
nated by sources of systematic error, notably the use of 
the quenched approximation, the lack of explicit extrap- 
olation to the continuum and possibly (we shall discuss 
this later) an overly small spatial volume. Nevertheless 
we still feel it is a valid exercise to compare the patterns 
seen in the data with the patterns present in the quark 



In addition, the meson source operator was placed a little fur- 
ther from the temporal Dirichlet wall, reducing the possibility of 
states propagating forward after reflection off the wall. 



model overlaps. The use of the quenched theory removes 
the possibility of multiquark states involving light quarks. 

As an illustrative example, we consider the T 1 chan- 
nel - at a finite lattice spacing correlators in this sym- 
metry channel will receive contributions from states that 
in the continuum have J PC = 1 ,3 ,4 .... With 
a set of twelve operators we extracted the mass spec- 
trum and overlap values shown in Table IIVI and Figure 
[5] By comparison with the quark model overlap forms 
we can propose a model-dependent interpretation of the 
spectrum that bares a strong resemblance to the conven- 
tional quark model picture. 

The ground state (of mass 3106(2) MeV) in this chan- 
nel has large overlaps with the unsmeared local opera- 
tors (7i,7o7i) indicating that it has a considerable wave- 
function at the origin and suggesting a dominantly 3 Si 
structure in line with what is expected for the J/tp. The 
first excited state (of mass 3764(18) MeV) also has large 
,707,) indicating 3 Si, and has a relatively sup- 
pressed coupling to the smeared local operators compared 
to the ground state - from the previous discussion we ex- 
pect this for a radial excitation and thus propose that 
this state is dominantly 2 3 Si. 

The second excited state (of mass 3846(12) MeV) has 
small overlaps onto all the operators used. A likely expla- 
nation of this is that we are getting overlap onto a 3 
state through lattice discretisation effects. In this case 
the overlaps would be suppressed by powers of the lat- 
tice spacing a relative to the overlaps onto the 1 states, 
which arc all that these operators have overlap with in 
the continuum (see the appendix of [ll|). Within the 
quark model a 3 ( 3 £>3) state is expected in this mass 
region and indeed in the T 2 ,A 2 ~ channels, where we 
include operators with continuum 3 overlap, there are 
candidate states whose Z values fit with belonging to a 
3 state (see [ll[ for details of the spin-assignment pro- 
cedure which depends only on continuum properties of 
the operators and not assumptions about models). 

The third excited state (of mass 3864(19) MeV) has 
small overlap onto the unsmeared local operators sug- 
gesting it is not dominantly a 3 Si state. It does however 
have a large overlap with the p x Dn operator which 
within the quark model has suppressed overlap onto 3 Si 
relative to 3 Di - this is due to the D operator transform- 
ing like Y 2 m . We propose that this state is mostly l 3 Z?i 
- its near degeneracy with the 3 state reflecting the 
small spin-orbit interactions in charmonium. It would 
appear that there is relatively little mixing between the 
nearby 2S, ID states, corresponding to a small tensor in- 
teraction. 

The fourth excited state (of mass 4283(77) MeV) has 
large overlap onto the unsmeared local operators suggest- 
ing that it may be dominantly 3 Si. We suggest that it is 
the 3 3 Si state. 

The fifth excited state (of mass 4400(60) MeV) has 
rather small overlaps onto all operators with the notable 
exception of the smeared ir x B, TT2 X B operators where the 
overlap is an order of magnitude larger than for the other 
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operator continuum allowed kinematic quark model 
name limit J P factor state 



7°7 5 



ijjijj 
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o+- 
o-+ 

o-+ 
1~ 



z 
z 
z 

z 

z e 
z e 

z € i 
z €< 



/(«) 



origin 
behaviour 



exotic 

x So 
3 Si 
a Di 
3 Si 
3 A 
3 Pi 





2^(1+4^ 



2V2 (l 
2^2 (l 



1 

4 m 7 



2 q" 
■ > mj. : 



2V2 (1 



2 

2 9 
\/3 rn q 



Pp(0) 


Ps(0) 
Ps(0) 
fls(O) 
Pfl(0) 
Ps(0) 
fl£(0) 
P'p(O) 
P-p(O) 



TABLE I: Local operators. Indicated are the J PC allowed at rest by Lorentz symmetry. The quark-model overlaps are given 
by Z = \/2M^n J tp(q)f(q) for unsmeared operators and Z^"' = \f2~M\pK j e~ CT 9 ^ 2 <p(q)f(q) for smeared operators. 
The final column indicates the leading behaviour in terms of derivatives of the spatial wavefunction at the origin. 
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kinematic 
factor 


quark model 
state 


no x V 






MZ 6* 


z s. 


«0(2) X V 


^7° a*v 


1-+ 


MZ €* 


exotic 


7T X V 




1+- 


MZ €' 


x Pi 


7T(2) X V 


^7°7 5 av 


l+- 


MZ €* 


l Px 


p x V 




0++ 


MZ5 ij 


3 Po 






1++ 


MZe ljk e k 


3 A 






2++ 


mz e tj 


3 P 2 
3 P 2 


P(2) X V 




0++ 


MZ& i] 


3 Po 






1++ 


Mze ijk e k 


3 P1 






2++ 


mz e lJ 


3 P 2 
3 P 2 


BixV 






MZ5 i] 


exotic 






f— 


MZe ijk e h 


3 Si 
3 i?i 






2~ 


mz e lj 


3 D 2 


61 x V 




0"+ 


MZ6 ij 


'So 






1-+ 


MZe ijh e k 


exotic 






2"+ 


mz e lJ 


'D 2 



f(q) 



origin 
behaviour 



2V2 g z 
3M m Q 

3M rriq 





4s/2 
y/3M' 
iV2 



iv3 fl (1- JL') 



V3M 



1 + 



1 

20m* 



5AJ m; 



3M y 



9 
4 /75 " 



75m 

4%/2 „ / , q 2 

\ 20m; 
J_ q 3 
5M m'i 



V3M 





2^2 q 2 
3M m 

3M mq 
2V2 q 2 
V5M rn q 
iV2 q 2 
3M m g 



8 q 2 
•/T&M m c 

TABLE II: As Table U for single derivative operators. 



R's(o) 



R' P (0) 
R' P (0) 
R' P (0) 
R' P (0) 
R' P (0) 
P'f(0) 
R'p(0) 
R' P (0) 
R' P (0) 
R'fW 


P's(0) 
P'd(O) 
i%(0) 
R's(0) 




states so far considered. The B-type operators have the 
characteristic feature that they are zero under the condi- 
tions of the quark model, corresponding as they do to the 
commutator of two covariant derivatives which vanishes if 
the gluonic field is neglected. With gluonic field included, 
the B operator is proportional to the field strength tensor 
which in a constituent gluon model would have at least 



one gluon creation/annihilation operator. Within an ex- 
tended quark model we would propose that considerable 
overlap onto B-type operators indicates some hybrid glu- 
onic nature to the state. In the flux-tube model, there 
is a degenerate set of exotics that includes a 1 state 
at around 4.2 GeV - within this model the hybrid 1 is 
a quark spin-singlet, which matches with the large over- 
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R'f(0) 

R'p(o) 


Rp(o) 



TABLE III: As Table U for two derivative operators. 



lap onto 7r x B which is dominantly spin-singlct. In the 
Coulomb gauge model [|[ there is a 1 hybrid state with 
a mass near 4.5 GeV which is also a quark spin-singlet. 

Although the masses do not agree terribly well with 
experimental candidates it would appear that the spec- 
trum we are observing is in line with the expectation of 
quark potential models with the intriguing addition of 
a state which may be hybrid in nature, perhaps being 
a candidate for the state being claimed in experiment 
at 4260 McV. We note here that as well as the masses 
appearing to be systematically high, we also overpredict 
the vector decay constants (proportional to Z^°\^i)) for 
excited states. Although we could not find any finite vol- 
ume effects in a limited study in [ll[ , we strongly suspect 
that this is the cause of these effects - if one numerically 
solves the Schrodinger equation with the Cornell poten- 
tial while forcing boundary conditions corresponding to 
finite volume and hence "squeezing" the wavefunction, 
one does observe systematically higher masses and larger 
wavefunctions at the origin high in the spectrum than 
the solutions in infinite volume. Future lattice work in 
larger volumes will clear this up. 

An example of a channel where the assignment is not so 
clear is , whose continuum content is h ,4 '", . . .. 



The large overlap, seen in Table [Vj onto the local opera- 
tors for all four states considered suggests they all have 
large 1 Sq components. Within the cc quark model this 
is the only possibility for pseudoscalar quantum num- 
bers. The near degeneracy of the second and third ex- 
cited states is not expected for radial excitations in the 
Cornell potential. The possibility of a lattice artifact 4 h 
( 1 (j ! 4 in the quark model, expected around 4.2 — 4.3 GeV 
with the Cornell potential) appears unlikely owing to the 
large overlaps with local operators. Within the Coulomb- 
gauge model of @ and the flux-tube model, there is a 
h hybrid expected degenerate with a 1 hybrid - we 
proposed such a state at 4400(60) in the T l analysis 
above. The second and third states suffer from large sta- 
tistical fluctuations - modulo this the large overlaps are 
with the smeared p x B operator and unsmcarcd local 
operators, suggesting hybrid and 1 Sq nature; the p x B 
operator is quark spin triplet as is expected of the hy- 
brid state. Within the large fluctuations there is room 
for some mixing of the conventional and hybrid states. 

Consideration of the channels (Ai,Ti,T2, E) ++ gives 
us an insight into the common properties of spin-orbit 
split multiplcts, e.g. 3 Pj. The spectrum presented here 
differs slightly from that published in owing to the 
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operator 


th [3106(2)] 1 


st [3746(18)] 2 


nd [3846(12)] 3 rd [3864(19)] 4 


th [4400(60)] 5 th [4283(77)] 


7 (»0 (10 -3) 


26.8(8) 


2.8(6) 


0.7(5) 


0.8(3) 


3.2(7) 


3.5(11) 


7 o 7 f m) (10- 3 ) 


26.2(8) 


3.3(6) 


1.1(5) 


1.3(4) 


2.8(8) 


3.5(11) 


a x V^ m) (HT 3 ) 


3.53(2) 


2.4(3) 


5.3(6) 


5.7(8) 


0.2(6) 


0.6(3) 


ai x V^ n) (10- 3 ) 


6.80(4) 


4.7(2) 


5.2(13) 


7.0(8) 


3.3(8) 


1.2(9) 


p x D^" 1 ' (10- 4 ) 


10.4(1) 


6(6) 


7(60) 


290(20) 


12(29) 


14(H) 


P 2 x D^ m) (HT 4 ) 


6.1(1) 


17(6) 


20(60) 


290(20) 


20(20) 


7(16) 


TT x B^r ) (1(T 5 ) 


30.6(5) 


4(5) 


3(7) 


3(6) 


130(10) 


25(25) 


7T2 X B^ m) (ID" 5 ) 


28.6(4) 


8(5) 


5(8) 


7(6) 


130(10) 


25(26) 


7i 


0.163(1) 


0.190(6) 


0.005(25) 


0.027(13) 


0.037(31) 


0.202(23) 


7o T< 


0.146(1) 


0.166(6) 


0.029(20) 


0.004(14) 


0.027(24) 


0.168(20) 


ai x Vti 


0.163(1) 


0.281(11) 


0.073(47) 


0.129(27) 


0.08(6) 


0.391(53) 


7T X BT1 (lO^ 3 ) 


8.89(9) 


8.5(5) 


1.4(10) 


1(1) 


8(3) 


8(3) 


7T 2 X Bti (10~ 3 ) 


9.09(9) 


8.0(5) 


0.7(10) 


0.5(8) 


8(3) 


8(3) 


assignment 


l 3 Si 


2 3 S'i 


3"- ? 


l 3 L>i 


hybrid ? 


3 3 Si 



TABLE IV: T" overlaps from numerical lattice QCD computation. Where an operator has the superscript (sm), the quark- 
fields have been Gaussian smeared - smearing parameters are given in Numbers in bold are characteristically large, 
indicating a particular quark model state assignment, while those in italics are both characteristically small and statistically 
compatible with zero within 3<r. The final row is a model-dependent state assignment described in the text. 



operator 


th [3027(2)] 1 


Bt [3714(27)] 


2 nd [4280(60)] 


3 rd [4500(150)] 


707^ m) (10- 3 ) 


24.2(8) 


7.0(9) 


0.41(14) 


7.6(8) 


7 (sm) (10 - 3) 


26.2(7) 


5.5(9) 


3.6(13) 


6.3(8) 


h x V^ m) (10- 3 ) 


9.46(5) 


5.9(5) 


7.4(21) 


6.3(8) 


p x B^ m) (10- 3 ) 


0.9(1) 


0.27(16) 


2.4(5) 


1.4(5) 


7075 


0.162(1) 


0.163(13) 


0.164(63) 


0.179(28) 


7s 


0.225(1) 


0.248(21) 


0.249(93) 


0.275(42) 


6l X V.41 


0.281(2) 


0.471(44) 


0.52(17) 


0.64(10) 


p x Mai (10 -2 ) 


3.00(3) 


2.6(2) 


0.8(17) 


3.5(6) 


assignment 


l 1 ^ 


2 1 5* 


hyb (+3^0) ? 


3 1 S (+ hyb) ? 



TABLE V: As [TV] but for the ,4"+ channel 



use of a somewhat larger operator basis causing an antic- 
ipated increase in the density of states produced through 
lattice discretisation artifacts. There is a clear ground 
state l 3 Pj=o,i,2 multiplet split by small spin-orbit forces 
and above that there is a dense spectrum of states whose 
overlap extractions are somewhat noisy; for reasons of 
space we do not show the numbers here, displaying only 
the mass spectrum in Figure [3] A plausible description 
would be that: 

• The A\ channel may house the 2 3 Po state at 
4080(340), although the noisy overlaps of this state 
are consistent with being a 4 ++ artifact; 

• The Ti channel houses the 2 3 P 1 state at 4119(59) 
and lattice artifact states at 4207(23), 4349(57) that 
might be the 1 3 F 3 (3++) and 1 3 F 4 (4++) states of 
the quark model; 

• The channel has four closely spaced levels 
around 4.1 GeV that have a possible interpretation 
as the 2 3 P 2 , l 3 F 2 , 1 3 F 3 , 1 3 F 4 states, only the last of 



which is a lattice artifact with the set of operators 
used; 

• The E ++ channel excited states near 4.1 GeV have 
large (but noisy) overlaps suggesting that they 
arc not lattice artifacts and are likely to be the 
2 3 P 2 and 1 3 F2 states. The reason for the non- 
appearance of a lattice artifact 4 ++ , seen in T\ and 
T2, is not known. 

We note that the overlaps onto the gluonic operator b\ x B 
arc uniform across all the states indicating within the 
quark model interpretation that there arc no dominantly 
hybrid states in this mass region - if our interpretation 
is correct this may favour the Coulomb-gauge model of 
heavy-quark hybrids over the simplest flux-tube model. 
In the Coulomb-gauge model [|| the lightest (0, 1,2) ++ 
hybrid states are heavier than the 1 hybrid state and 
are quark spin-triplets, while in the simplest flux-tube 
model there is a spin-singlet 1 ++ hybrid state degener- 
ate with the 1 hybrid state, which we have proposed 
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4600 



4400 - 



4200 - 



4000 " 



> 

a> 

^ 3800 

in 
in 

j= 3600 h 



3400 - 



3200 " 



3000 " 



hyb 
4>S 



— u, 

hyb?f 




i| 

4400 


3 3 S, 









3 3 S, = = = = ! 



i)<(4040) 



1 3 n, l I 3864±1 9 

1 3 D,L-_-_-_-_-J i]i(3770) 



2 3 V 







ij)(3686) 



1 3 S,' 



1 J S, 



JA|) 



FIG. 2: 1 __ spectrum. Black lines are experimental states 
at various levels of confirmation, blue bars are lattice states 
with quark model assignments as described in the text. Also 
shown are quark model predictions and the hybrid prediction 
of the Coulomb gauge model. 



is located at a mass near 4.4 GeV, and since the b\ x B 
operator has a quark spin-singlet component it seems rea- 
sonable to anticipate overlap with this state. 

Finally we will consider applying our methods to chan- 
nels that in continuum house J PC exotic states not 
present in the simple cc quark model. The most straight- 
forward choice is which only houses even spins and 
hence is always exotic with this PC. The simplest op- 
erator transforming in this way is ■07° tp - we find that 
unsmeared correlators containing this operator are con- 
sistent with zero, which matches with the fact that this 
operator contains no gluonic field to generate a hybrid 
state. By smearing this operator using a function of the 
gaugc-covariant laplacian we introduce the gluonic field 
and have the possibility of overlap with a hybrid state. 
Using this operator and smeared and unsmeared a\ x B 
operators we extract a state at 4465(65) MeV which has 
clear overlaps on to all three operators. The first excited 
state is much higher in mass, at 5570(270) MeV. 

A less trivial channel is T 1 _+ which houses exotic 
1 h , 3 h but also non-exotic 4 h which is realised in the 
quark model as 1 C±. In [ll[ it proved to be not possible to 
decisively state whether the ground state in this channel 
was indeed the exotic 1 h state or a 4 ^ nonexotic, here 
we will consider this again using our model-dependent 
overlap comparison. In Table PVTl we present the extracted 
overlaps using an eight dimensional basis of operators. 
For the ground state we see rather large overlaps with 
the quark-model forbidden operators p x B, p 2 x B sug- 
gesting that it may well be the exotic hybrid 1 h state. 
The first excited state has overlaps consistent with zero 



for all operators except the smeared b\ x V and possibly 
the a (2) x V. This may well be a signal for a lattice 
artifact 4 h state, whose spin-singlet ^Ga) nature could 
appear through overlap with b\ x V which has a dominant 
spin-singlet piece in its non-relativistic reduction. 

2 H exotics appear in (T 2 ,Ey . The T 2 receives 
contributions also from non-exotic 3 H and indeed the 
ground state in that channel is identified as such. The 
first excited state in T 2 matches with the ground state 
in and in both cases large overlaps with the a\ x B 
operator are seen strongly suggesting that this is the 2 +_ 
exotic at a mass of 4620(60). 



V. FINITE MOMENTUM 

The non-relativistic quark-model states as constructed 
do not transform covariantly under boosts, but do trans- 
form properly under rotations in three-dimensions. This 
can lead to overlap on to more states than are allowed 
by Lorentz symmetry. For example, consider the opera- 
tor tpj^d^ip - insisting upon Poincare invariance one has 
only the following overlaps 

(0|^9^|0 ++ (p,r)) = Z^ + Z'pV 
(0|^7 AI ^|l++(p,r)) = Ze^p p e a {p,r) 
(0|VV I 1 9 1 >|2+ + (p») = Ze"»), 

so that if, as we do in the lattice calculation, one considers 
only the spatial derivatives, one has overlaps 

(0\4>Yd j tp\0 ++ {p,r)) = Z6" ■ Z'p'p 1 

(0|V5 7 l a j V|l ++ (p,r)) = Z^ k (p k €°(p,r) - E e fe (p») 

(0|Vi7 l a j V|2 ++ (p,r)) = Ze i3 (p». 

But note that this is not the most general set allowed by 
three-dimensional rotations, parity and charge conjuga- 
tion, giving up on boost invariance we also are allowed 
overlaps 

Ze^ k p k 



(0\^d j i}\Q- + (p,r)) 
(0|^9^|2- + (p,r)) 



(p, r)p l 



Within the quark model state construction, we can ex- 
plicitly compute these overlaps at finite momentum find- 
ing 



Z( 3 P ) 
Z'( 3 P ) 



ZfPx) 
Z( 3 P 2 ) 
ZCSo) 

zCd 2 ) 




q 2 dq 2 
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FIG. 3: PC = ++ spectrum by lattice irrep and by continuum assignment. 



operator 


th [4305(40)] 


1 st [4645(86)] 2 nd 


[4689(138)] 3 rd 


[5580(160)] 


o 0( a) x V^ m) (10- 3 ) 


2.5(2) 


2.0(6) 


2.0(7) 


0.8(4) 


bi x V^ m) (10- 3 ) 


2.2(2) 


2.9(4) 


1.7(8) 


1.4(7) 


p x B^ m) (10~ 3 ) 


2.88(5) 


0.2(2) 


0.8(5) 


0(0.3) 


p 2 x By" 1 ' (1(T 3 ) 


2.84(5) 


0.0(2) 


0.8(5) 


0(0.3) 


a (2) x Vti (10~ 3 ) 


1.8(1) 


0.2(3) 


1.2(5) 


1.5(7) 


bi x Vti (10 -3 ) 


1.7(1) 


0.2(4) 


1.0(5) 


0.8(11) 


p x B T i (10~ 3 ) 


3.1(2) 


0.4(7) 


2.6(9) 


3.6(20) 


pa x B T i (10^ 3 ) 


3.0(2) 


0.2(6) 


2.5(9) 


3.5(20) 


assignment 


1"+ hyb ? 


4-+( 1 G 4 ) ? 


? 


? 



TABLE VI: As [TV] but for the Tf+ channel.. 



so that, as anticipated, as well as the overlaps allowed 
by Lorentz symmetry there are also disallowed over- 
laps. This is an inherent weakness of the non-relativistic 
model that can only be remedied by constructing a fully 
Poincare covariant bound-state scheme which poses a sig- 
nificant challenge to modellers. 



VI. SUMMARY 

We have presented a simple framework for compari- 
son of lattice QCD spectroscopy and the quark model. 
It relies upon a non-relativistic reduction and as such is 
suitable in the heavy-quark sector. We have compared 
with recent lattice QCD data and presented a model- 
dependent description of the data which agrees in struc- 
ture with the predictions of the Cornell potential quark 
model, but goes beyond that model in providing predic- 
tions for exotics and crypto-exotic hybrid mesons. The 
particular lattice data used is dominated by sources of 
systematic error, notably we suspect that the small vol- 
ume used is "squeezing" the wavefunctions of higher ex- 
cited states. An alternative use of this method, when 
applied to more realistic lattice QCD data, would be to 
allow a quark model to be "tuned" to QCD, through se- 
lection of interactions and parameters. 



Of particular interest is the extension beyond quark- 
model states, where the gluonic field plays a manifest 
role. If the assignments of hybrid nature made in this 
analysis are correct we have the following (incomplete) 
hybrid spectrum: 

• a non-exotic pseudoscalar state (0 h ) at 4280(60) 
MeV which may have a degree of mixing with a 
nearby conventional cc state; 

• an exotic 1"+ state at 4305(50) MeV; 

• a non-exotic vector state (1 ) around 4400(60) 
MeV where mixing with conventional states is not 
apparent; 

• an exotic 0+~ state at 4465(65) MeV; 

• an exotic 2+~ state at 4620(60) MeV; 

• no non-exotic hybrids in (0, 1, 2) ++ channels below 
about 4.5 GeV 

An interesting extension to the work done in this paper 
might be to apply a similar technique using a model with 
explicit gluonic degrees-of-freedom such as the flux-tube 
model or the Coulomb-gauge model. 

It is not clear if this method will have utility for lighter 
quarks, where the quasiparticlc quark- like degrees-of- 
freedom in the quark model ("constituent quarks") are 
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not the same quark degrees-of-frccdom that appear in the 
QCD Lagrangian ("current quarks"). 
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